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We consider the effect of a six-fermion interaction of the 't Hooft form in the quark-quark channel 
on the ground state of matter at finite density. The coupling constant for this new term is varied 
within the limits suggested by naturalness. The flavor-mixing effects of the additional term desta- 
bilize the color-flavor-locked (CFL) and, to a lesser extent, the two-flavor color superconducting 
(2SC) phases of quark matter, especially for positive values of the coupling. For some values of 
the coupling, the critical density for CFL phase is nearly larger than the maximum density in the 
neutron stars. We comment on the implications for neutron star evolution. 
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PACS numbers: 12.38.-t, 12.39.Fe, 26.60.+C, 97.60.Jd 



HISTORY AND MOTIVATION 

The study of dense (baryon chemical potential ~ 
1.5 GeV) matter has been recently revolutionized by 
the observation that dense quark matter exhibits color- 
superconductivity and that the gaps may be of order 100 
McV Q, Q) El • Gaps of this magnitude are large enough 
to have significant implications for neutron star struc- 
ture 0|, proto-neutron star evolution and neutron 
star cooling 0. 

Since directly utilizing QCD at the relevant densities 
(Mbaryon ~ 1 GeV) is so far impossible, the use of effective 
theories like the Nambu-Jona-Lasinio (NJL) model is 
common in the study of dense quark matter. In the NJL 
model, the high-energy degrees of freedom (the gluons) 
are integrated out and we restrict ourselves to workin 
at energy scales less than the momentum cutoff A 0, 



-Ceff = 



L^i Adim(0 ra )-4 



(1) 



where O n are operators, dim(C„) is the dimension of the 
operator, and c„ are dimensionless coupling constants. 
Because it is impossible to create a model with the same 
symmetries as QCD with four-fermion operators alone 
the so-called 't Hooft term 



C 



-"kinetic 



H~ £four— fermion H~ tHooft 



tHooft ~ det f [q(l + To)q] +det/ [q(l -75)9] (2) 

is added where q is a quark spinor and det t is a determi- 
nant in flavor space. This term, like QCD, respects chiral 
symmetry but breaks the U(1)a symmetry. The use of 
this Lagrangian is the standard approach which has been 
used to describe dense matter. 

Unfortunately, when employed to study quark super- 
conductivity, this standard approach does not employ a 
manifestly consistent truncation scheme; the quark-anti- 
quark interaction is treated at the six-fermion level, but 
the quark-quark interaction is only treated at the four- 
fermion level. There is no reason to rule out a term of 
the form (omitting color for clarity) 

£ C S6 ~ K mQ e ijm eun (<?»75«f )(qf 75%) (<7m?n) (3) 



which has the same symmetries as C> tHooft- This term 
may have a significant impact on the nature of dense 
matter 0, H El . While the effect of the dynamically 
generated quark masses on the superconducting gaps has 
been studied 0] , Eq. |31 implies a modification to the 
quark masses due to the presence of the gap. 

In this article, we study the effect of £cS6 on the quark 
masses and on the phase structure of dense matter. We 
show that sufficiently positive values of ATdiq increase the 
quark masses and thus favor less-gapped phases, while 
sufficiently negative values of i^DiQ split the magnitude 
of the up-down and light-strange gaps. The phase struc- 
ture of dense matter thus depends critically on the sign 
and magnitude of this unknown parameter. 



THE MODEL LAGRANGIAN 

The Lagrangian is [13, 13 



C 



a=0 



-K [det f q (l - i-f) q + det/g (l 

X (qtail5qff3)(qk-y i l5q^)(q m sqn n ) 



• i 7 5 ) q] 



(4) 



where flavor is represented by Latin indices, color is 
represented by Greek indices, and the charge conju- 
gate Dirac spinors are defined by (Cj^C — 7 pT and 
C T = -C) 

q C = Cf and f = q T C . (5) 

The four-fermion coupling in the quark-anti-quark chan- 
nel is denoted G, the four-fermion coupling in the quark- 
quark channel is denoted Gdiq, and mo.ij is the constant 
current quark mass matrix which is diagonal in flavor. 
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We utilize the ansatze 

gi<?2<Z3<?4 (?ig2> q3Q4 + qiq2 (9394) - (9192) (9394} 

9l 9293949596 9l92 (9394) (95%) + (?1?2> (fe^) 95% 

+ (9192} 9394 (qsqe) 

-2 (gx^) (9394) (f5?6> (6) 

to obtain the mean-field approximation |l5[ . This pro- 
cedure retains only the lowest order terms in the 1/N C 
expansion ^| . Color neutrality is ensured using the pro- 
cedure from Ref. [l^ . In the mean- field approximation, 
the Lagrangian is only quadratic in the fermion fields, 
and the thermodynamical potential can be obtained from 
the inverse propagator in the standard way (it| . The mo- 
mentum integrals in the gap equations are divergent, and 
are regulated by a three-momentum cutoff denoted by A. 
The inverse propagator is numerically diagonalized for 
each abscissa of the momentum integration to obtain the 
thermodynamical potential. 

For simplicity, we sometimes use the notation ~ 
e ijk ^qC^5q^^ SQ gaps are denoted with the flavor 
of quark that is not involved in the pairing e.g. A uc i 
is denoted by A s . Other than {qiqi) and A,, we as- 
sume that all other condensates vanish. This includes the 
pseudoscalar condensates which which are likely present 
in dense matter and naturally accompany the Goldstone 
bosons [l8|. This (non-trivial) complication will be left 
to later work. 

The effects of £cs6 on the thermodynamical poten- 
tial can be summarized in three modifications from the 
standard approach where ATdiq = 0- These changes are 
that the values of the gap in the inverse propagator are 
modified 



A, A,(l + -^-(g^)) . 

-iV c (_zDIQ 



(7) 



a new effective mass term (which includes contributions 
which are not diagonal in flavor) appears 



ATdiq _ a a 

—2 9i A i^j1j 



AG 



(8) 



DIQ 



which modifies the dynamical mass 



m , = mi >0 - 4G (qq) t + K\e ijk \ (qq) j (qq) k + ^£3_A?(9) 

4G DIQ 



and that there is a new contribution to the vacuum en- 
ergy 



if DIQ 

^diq = — y . A i w<?». 

ZCr DIQ 



(10) 



Note that Eq. [SJmeans that the quark masses are density- 
dependent if ATdiq 7^ even when the chiral conden- 
sates (qq) vanish. As has been suggested [l^, this term 



generates a dynamical quark mass entirely distinct from 
the typical mechanism of spontaneous chiral symmetry 
breaking. The quarks obtain a dynamical mass even 
when the quark condensate is taken to be zero. 

One may use Fierz transformations to calculate the 
coefficient ATdiq from the quark-anti-quark form of the 
't Hooft interaction. One can view this in the follow- 
ing way: For each prospective new term, e.g. the term 
M7 5 d c d c 7 5 MSS, there are six Fierz transformations (for six 
fermions this is a 35 x 35 matrix instead of the usual 5 
x 5 matrix for four fermions) that give this term when 
applied to the six sets (in flavor space) of terms in the 't 
Hooft interaction. One such transformation is 



qiqjqkqeq m q n + qci q 3 qkj qs.q m q n 
-+ 
l 



+qiqjqki b qtq m i 5 q n + qa 5 qjqkq£q m ^ 5 qv 



= 2 U»7 5 ?k ?f 7 5 «e<?mfti + qiqkqfi qiqr,ii b qn 
+qtqkJ 5 qfqeq-mi 5 qn + q%q k qf qtq m qn) 



■- (q l -f 5 o-^qkqjJ 5 o- fJ ,uqeq m qn 

-qa 5 <r^ v qkqf <r»v<itq~mi 5 qn) 



(ii) 



When these transformations are combined to give the 
coefficient i^DiQ, the result is zero (see the Appendix). 
Although terms with different Dirac structure do survive, 
we do not include these terms since we do not expect the 
corresponding condensates, e.g. (qa^q) to be non-zero. 
This procedure is not the only possible approach for de- 
riving the mean-field Lagrangian (one could also enumer- 
ate all possible Wick contractions). Because alternate 
approaches and/or using terms of higher order in 1/N C 
may modify this result, we cannot conclude necessarily 
that if diq must be zero. Also, terms like Eq. [3] with a 
Dirac structure {qq)(qo-^q)(qo' fJ ' l 'q) H3 and their corre- 
sponding terms in the quark-quark channel may play a 
role. We leave these considerations to future work. 

Our Lagrangian is free to contain any terms which fol- 
low the symmetries of the underlying theory. We expect 
that the coefficient of this term will be "natural" . When 
the coefficients are expressed in terms of the underlying 
length scales (in our case, the momentum cutoff A), the 
coefficients should all be of similar magnitude. We al- 
low the coefficient ATdiq to vary, between the values —AT 
and AT, which we view to a modest variation as suggested 
by the constraints of naturalness. A larger variation in 
i^DiQ is not necessarily excluded. We use the values of A 
and the current quark masses from Ref. [HI > where they 
are fixed by matching the pion, kaon, and rf masses in 
vacuum as well as the pion decay constant. We choose 
to fix GdiqA 2 = 1.61 to be large enough so that the 
maximum value of the gaps (when including the quark 
dynamical mass) as a function of density when ATdiq = 
is about 80 MeV, close to the value of about 100 MeV 
predicted by calculations in perturbative QCD. If the 
dynamically-generated quark mass is assumed to be zero 
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and the mass-gap equations are ignored, then the maxi- 
mum value of the gap predicted by this model is about 
120 MeV. We leave Gdiq fixed when varying ifoiQ- One 
could also allow the coefficient, Gdiq as a function to 
vary as a function of -Kdiq by instead ensuring that the 
maximum value of one of the three gaps at high densities 
is constant. We have checked that this alternative does 
not change our conclusions significantly. 



RESULTS 

Obtaining analytical results is difficult, due to the fla- 
vor mixing mass terms from Eq. [S] which make it diffi- 
cult to directly reduce the inverse propagator (a 36 x 36 
matrix) into a block-diagonal form. It is possible, with 
High-Density Effective Theory > to simplify the Dirac 
structure, but this would likely result in a 9 x 9 inverse 
propagator which is also difficult. It is also possible to 
restore the usual form of the propagator encountered in 
studies where -Kdiq = by ignoring the terms in Eq. 
|H1 where i ^ j. In this case, the inverse propagator is 
worked out in detail in Ref. [23j. 

It is possible to see qualitatively, what the effect of 
adding a term with -Kdiq ^ might be from Eq. 
above. When -Kdiq > 0, we expect the gaps decrease as 
the quark condensate increases, and thus the gap should 
decrease with increasing mass. However, from Eq. [3] we 
expect the opposite and we find that it is this effect that 
dominates the description of the strange quark mass and 
A u d- Further complicating the analysis, Eq. 03 indicates 
that an increase in A us and Ad s will tend to split the 
mass of the up and down quark, thus possibly weakening 

Aud- 

We study charge- and color-neutral, beta-equilibrated, 
bulk matter at fixed baryon density and a fixed temper- 
ature of 10 MeV. We operate at a small but finite tem- 
perature in order to alleviate the numerical difficulties 
of discontinuities in the momentum integral present in 
the thermodynamical potential. The zero-temperature 
results will not deviate significantly from our results. We 
include non-interacting electrons, but we do not include 
neutrinos. The addition of neutrinos would further split 
the approximate flavor symmetry between the up and 
down quarks. Our results will faithfully describe mat- 
ter in the center of a neutron star containing quarks a 
minute or later after formation [2"! l2fl|. 

We note that the effect of -Kdiq is small when the 
quark condensates (qq) are taken to be zero. When this 
is assumed to be true, then the modifications from Eqs. 
and[H]have no effect, and the gaps are nearly independent 
of -Kdiq . We remove this assumption and solve the mass 
gap equations for the quark condensates in the following. 

Figure ^ presents the masses and gaps in the CFL 
phase at fixed density and temperature as a function 
of Kdiq. Both the quark masses and A u d increase as 
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FIG. 1: Quark masses (dotted lines for u and d, and dashed- 
dotted line for s) and gaps (solid line for A u d and dashed line 
for the other two gaps) as a function of -Kdiq in the CFL 
phase at ub = 1.6 fm -3 . 



Kdiq increases. The effect from Eq. [5] causes the A us 
and A^ gap to decrease when A u d increases. At suffi- 
ciently large values of the coupling, the strong increase of 
the strange quark mass destabilizes the CFL phase. For 
Kdiq > 0.4, the gap equations have no solution. If the 
coupling was verified by some other means to be larger 
than this critical value, then the CFL phase could not be 
present at this density. As Kjuq/K — > —1, the effects 
on the masses and gaps tend to be less extreme. The 
most significant effect is the increasing split between the 
values of the light-strange gaps, A us and Ad s , and the 
light-quark gap A u d- One might expect the dependence 
of the gaps on -Kdiq would change the phase structure of 
matter by shifting the energy density. However, we find 
that this is not the case and that the energy density is 
relatively constant as a function of -Kdiq- Note also that 
the strange quark mass can change by as much as 50% 
for different values of -Kdiq- 

Also plotted in Figgis the parameter to^/(/iA) where 
/i and A in this context are computed by averaging the 
quark chemical potentials and gaps over all flavors and 
colors. This parameter has been demonstrated to be the 
relevant dimcnsionlcss quantity which dictates the phase 
content of quark matter at high density r 26|. Values of 
mf/(/iA) larger than 2 suggest a transition to a gap- 
less CFL phase jifl ]. while values larger than 4 suggest 
a transition to the 2SC phase. The presence of the gap- 
less 2SC phase [2^] is also probable when .Kdiq becomes 
positive, since the ud gap is becoming stronger and the 



4 



light-strange gaps are weakening. The value of m 2 s / (//A) 
is also important in dictating the number and type of 
Goldstone bosons present in the CFL phase. This pa- 
rameter is quite flat for small variations in Kbiq, but 
increases or decreases strongly when the absolute magni- 
tude of i^DiQ is sufficiently large. When Kbiq/K ~ 0.4, 
the value of mg/(/iA) is nearly 2, indicating the disap- 
pearance of the CFL phase. 




FIG. 2: Strange quark mass (dashed-dotted line) and the gap 
(solid line) as a function of -Kdiq in the 2SC phase at ng = 1.6 
im" 3 . 

The results for the 2SC phase, where A us = Ads = 0, 
at the same density and temperature are shown in Figure 
12 The strange quark mass and gap increase strongly for 
increasing Kbiq as in the CFL phase, leading to a critical 
value above which the gap equations have no solution. At 
this density, the 2SC phase is also not present in matter 
for Kbiq > 0.4. The similarity of this critical value of the 
coupling to the CFL phase in Fig. ^ is a result of the fact 
that the light-strange gaps and light-quark masses are 
comparatively small and thus do not strongly perturb m s 
and A u d- The parameter mf/^A) is slightly modified 
from the CFL case and indicates that the 2SC phase 
is also likely to be unstable for large negative values of 
Kbiq/K. 

For comparison, Fig. presents the quark masses and 
gaps in the CFL phase at a lower density, 1.2 fm -3 . The 
results are not much different from Fig. ^ The quark 
gaps have decreased slightly and the quark masses are 
slightly larger. The major distinction is that the criti- 
cal value of i^DiQ above which the CFL phase does not 
appear has been lowered from 0.4 to less than 0.3. 

The implication of the shift in the critical value of 
-Kdiq is that the critical density for the onset of the 
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FIG. 3: The quark masses and gaps in the CFL phase at 
us = 1.2 far 3 . See Fig.0 
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FIG. 4: The gaps as a function of baryon density in the CFL 
phase for two values of the coupling A'diq. 



CFL phase is drastically affected by a positive value of 
Kviq- This is demonstrated in Fig.0] where the smallest 
quark gap vanishes at about 1.0 fm~ 3 when Kbiq = 0, 
and at about 1.42 fm -3 when Kbiq/K — 0.3. On the 
other hand, because the gaps are not as strongly mod- 
ified when i^DiQ is negative, the critical density when 
Kbiq/K = —0.5 is almost unchanged, moving down only 
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to 0.96 fm -3 . In this figure, the Ad s gap does not vanish 
completely to zero because the solution of the gap equa- 
tions becomes difficult when the gaps are small. Note 
again that an increase the parameter m 2 s / {p,A) indicates, 
to some extent, the disappearance of the CFL phase as 
the density decreases. 

DISCUSSION - NEUTRON STARS WITH 

KniQ + 

These results may have several implications for neutron 
star structure and evolution. 

Neutron-star masses and radii: We have solved the 
Tolman-Oppenheimer-Volkov equations for neutron star 
structure under the assumption that there is no mixed 
phase (i.e. the surface tension is very large so that mixed 
phases are suppressed). The results for Kbiq = and 
Kdiq/K — 0.3 are given in Fig. |5] We find that neutron 
star masses and radii are not very sensitive to i^DiQ for 
the model that we have chosen. As mentioned above, a 
positive value of -Kdiq tends to increase the critical den- 
sity for the appearance of the CFL phase, thus stiffening 
the equation of state and slightly increasing the maxi- 
mum mass from 1.83 to 1.9 Mq. The small magnitude of 
this effect is partially due to the fact that, for Kdiq = 0, 
the critical density for the appearance of quark matter is 
1.0 fm~ 3 , while the central density of the maximum mass 
neutron star is only 1.45 fm~ 3 . There is not much quark 
matter to begin with, so therefore the effect of -Kdiq is 
limited. In regards to the phase content of matter in the 
neutron star, the larger value of -Kdiq nearly pushes the 
CFL phase out of the neutron star entirely, and the cen- 
ter of the neutron star is dominated by the 2SC phase 
instead. These effects will be larger if the diquark cou- 
pling is increased and may be modified by the presence 
of a mixed phase. Recent calculations of the surface ten- 
sion suggest that it is small, and thus a mixed phase may 
be present (2^. It would be interesting to examine the 
effect of -Kdiq on this surface tension. 

Neutron-star cooling: Neutron stars are sensitive to 
the difference between the CFL and 2SC phases since 
the former is likely to contain no 

(or very few) elec- 
trons. Neutron star cooling is affected by the presence 
or absence of electrons because the specific heat of mat- 
ter is dominated by the electrons when they are present. 
In addition, the specific heat contribution from the light 
quarks is proportional to exp(— Ag s /T), where £ — u or d 
which is exponentially small in the CFL phase and of 
order unity in the 2SC phase. Also, the splitting of the 
gaps at Kbiq/K = — 1, will enforce two critical temper- 
atures for quark matter in the CFL phase. The so-called 
"gapless" phases are dependent upon the strange quark 
mass which is strongly affected by if diq , especially when 
it is greater than zero. These gapless phases, because of 
the nature of the quark dispersion relations, have unusual 
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FIG. 5: The mass versus radius for neutron stars with two 
values of the coupling -Kdiq. 



transport properties that also have implications for neu- 
tron star cooling [2H El l3l| . 

Proto-neutron star evolution: Ref. [f| pointed out that 
the proto-neutron star neutrino signal may be increased 
noticeably by the enhanced cooling that is present when 
the core temperature falls below the critical temperature. 
The presence of -Kdiq implies that this enhanced cooling 
may occur in (at least) two separate stages, as the critical 
temperatures corresponding to the A u d and the light- 
strange gaps are surpassed. 

To the extent that the presence of a large (~ 100 MeV) 
gap affects observations of neutron stars, the presence of 
a color-superconducting six-fcrmion interaction also has 
an impact on neutron star observations. It would be 
interesting to compare these results with those from gap- 
less CFL phases. Because of the uncertainty in the value 
of the -Kdiq coupling (in addition to the other uncer- 
tainties already present), it will be difficult to settle the 
questions of the nature of dense matter and the question 
of the properties of neutron stars containing deconfined 
quark matter until theoretical progress in QCD or as- 
trophysical observations can more completely settle the 
issue. 

The author would like to thank T. Bhattacharaya, T. 
Biirvenich, M. Prakash, S. Reddy, G. Rupak, and T. 
Schafer for discussions. This work has been supported by 
the DOE under grant number DOE/W-7405-ENG-36. 
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Appendix A - Fierz Transformations 

Traditionally, the Fierz transformation is denned as 
the matrix, C, which obeys the relation 

T [^iTa^j^uTkfl)^ = ^2 C a ,bi>iTb,ieipe'ipkF b kj ipj (12) 



where I\ 

|"l,7 5 ^ 'V 5 'T/^ 



[1,7 5 >7/*>7 5 7m><v] and pi e 
,7 f ",7"7 f ",(T f "* / ] for i = a, 6. Because the Fierz 
transformation is nothing other than a set of equalities, 
the two four-fermion interactions on both sides are 
necessarily equivalent. However, in the mean- field 
approximation, these two forms lead to different thermo- 
dynamic potentials. For simplicity, the transformation 
above will be denoted J-(ijkl — > itkj). 

One may also perform a four-fermion Fierz transforma- 
tion in the quark-quark channel, namely, J-(ijk£ — ► ikji) 
(see the review in Ref. 32] ). This transformation op- 
erates over a different "basis" of combinations of Dirac 
matrices 



(13) 



where e [75, 1, 7 5 7m> 7h, °>] rV 6 

[7 s , 1, 7 5 7 At , 7 M , cr A1I/ ] (remember that ip c j 5 tp is a 
Lorentz scalar). We can simplify the notation for 
the basis by using the notation of a direct product: 
7 5 <g> 7 5 , 1 <8> 1, 7 5 7 M ® 7 5 7 P , 7 M ® 7 M , o> ® cr' 1,y , or more 
simply, SS, PP, VV, AA, TT. 

In a similar notation the 35-element basis for the six- 
fermion Fierz transformations is 

SSS, SPP, PSP, PPS, SVV, VSV, VVS, 
SAA, ASA, AAS, STT, TST, TTS, 
PVA, PAV, VPA, VAP, APV, AVP, 
TVV, VTV, VVT, TAA, ATA, A AT, TTT, 
VAQ, VQA, AVQ, AQV, QVA, QAV, 

PTQ,TPQ, QTP (14) 

where Q denotes a "pseudo-tensor" combination, 7 5 cr MI/ . 
(The Q terms are used as an alternative to the formula- 
tion in terms of objects of the form E K \ ilv a llv 33J.) All 
other combinations can be written as a linear combina- 
tion of these 35 basis elements. In order to distinguish 4- 
and 6-fermion transformations, we will use a subscript, 
i.e. ^{ijkimn — > iikjmn) is really a four-fermion trans- 
formation since the fields with indices m and n are not 
participating in the transformation. 

The Dirac scalar terms in the 't Hooft interaction are 



uuddss + usdusd + uddssu 
—udduss — usddsu — uudssd 



plus the corresponding terms created by adding an even 
number of 7 s matrices. 

As a demonstration, we examine the coefficient of the 
term u r y 5 d c u c j 5 dss. The various contributions to this 
coefficient are 

+J-i{ijklmn — > ikjlmn) [uuddss + uudj 5 dsj 5 s 

+U7 5 udds'f 5 s + uj 5 udj 5 dss\ 
+J 7 Q(ijklmn — > iklnmj) [usdusd + usdj 5 usj 5 d 

+uj 5 sdus'f 5 d + vrf sdy usd] 
+J-§{ijklmn — > iknjml) [uddssu + uddj ss'y u 

+uj 5 ddss'f 5 u + uj 5 ddj 5 ssu] 
—Tt±{ijklmn — > ikljmn) [udduss + udd^us'j 5 's 

+uj 5 ddus'f 5 s + vrf ddry uss] 
—J-(i(ijklmn — > iknlmj) [usddsu + usdj 5 dsj 5 u 

+uj 5 sdds / -f 5 u + uj 5 sdj 5 dsu] 
—J-(i(ijklmn — > ikjnml) [uudssd + uudj 5 ss"f 5 d 

+uj 5 udss'f 5 d + uj 5 udj 5 ssd] 

Note that result of the first of these six transformations 
is given as Eq. in the text. The coefficients of the 
desired term, u"f 5 d u c ^ 5 dss, from each of the six trans- 
formations (together with an appropriate factor of -1 for 
odd fermionic permutations) are 



111 

"2> 4> 4 



(16) 



and the sum is zero. 

Because the coefficient is zero, we need not consider 
the Fierz transformations in the SU(3) (color or flavor) 
spaces. However, since the result for six-fermion trans- 
formations in SU(3) is not present in the literature, we 
give the 15-element basis for computing the transforma- 
tions (this enlarged basis from that presented in Ref. j^il 
is necessary to perform the transformations in the quark- 
quark channel) 



ibc 



1 <g) 1 

1® A a 

A a <g> \ b < 



;1, 

> A a 
A c 



A a ® A a <g> 1, 

A a <8> 1 <g> A a , 
ifabcX a ® A b ® A c , 
X(A) a ® X(A) a ® 1, 1 ® X(A) a (g> X(A) a , 
X(A) a ® 1 ® X(A) a , 
\{A) a \{S) b ®\(A) a ® X(S) b , 
X(A) a X(S)b <8 X(S) b <& X(A) a , 
X(S) b ® X(A) a X(S)b <8> X(A) a , 
X(A) a ® X(A) a X(S)b ® A(S') b , 
A(A) a ® A(S) 6 ® A(A) a A(S')b, 
A(S') b ® A(A) a ® A(A) a A(S')b 



(17) 



(15) 



where the occurrence of A indicates an implicit sum over 
all 8 SU(3) matrices, X(A) restricts the sum to only the 
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three anti-symmetric A matrices, and the implicit sum 
for X(S) is over the six symmetric A matrices. The full 
results are available from the author. 
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